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Variation of the Amati Relation with the Cosmological 
Redshift: a Selection Effect or an Evolution Effect? 
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ABSTRACT 

Because of the limit in the number of gamma-ray bursts (GRBs) with available red- 
shifts and spectra, all current investigations on the correlation among GRB variables 
use burst samples with redshifts that span a very large range. The evolution and se- 
lection effects have thus been ignored, which might have important influence on the 
results. In this Letter, we divide the 48 long-duration GRBs in Amati (2006, 2007) 
into four groups with redshift from low to high, each group contains 12 GRBs. Then 
we fit each group with the Amati relation log-Ei SO = a + 61ogi?p Ca ] < , and check if the 
parameters a and b evolve with the GRB redshift. We find that a and b vary with the 
mean redshift of the GRBs in each group systematically and significantly. Monte-Carlo 
simulations show that there is only ~ 4 percent of chance that the variation is caused 
by the selection effect arising from the fluence limit. Hence, our results may indicate 
that GRBs evolve strongly with the cosmological redshift. 
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m 1 1 INTRODUCTION 



A remarkable progress in the observation of gamma-ray 
bursts (GRBs) has been the identification of several very 
good correlations among the GRB observables (see Schaefer 
2007 for a review). Based on several of those correlations, 
some people have eagerly proposed to use GRBs as stan- 
dard candles to probe the cosmological Hubble diagram to 
very high redshift (Schaefer 2003; Dai, Liang & Xu 2004; 
Ghirlanda et al. 2004b; Lamb et al. 2005; Firmani et al. 
2006; Schaefer 2007, and references therein). Enlightening 
comments and criticism on GRBs as standard candles can 
be found in Bloom, Frail & Kulkarni (2003) and Friedman 
& Bloom (2005). 

All the GRB correlations have been obtained by fitting 
a hybrid GRB sample without discriminating the redshift. 
Indeed, the redshift in the sample usually spans a very large 
range: from z ~ 0.1 up to z ~ 6. This is of course caused 
by the fact that we do not have an enough number of GRBs 
with measured redshifts limited in a small range. Then, in- 
evitably, the effect of the GRB evolution with the redshift, 
and the selection effects, have been ignored. This raises an 
important question about whether the relations that peo- 
ple have found reflect the true physics of GRBs or they are 
just superficial. [See Band & Preece (2005) for a nice discus- 
sion on the selection effect and the correlation between the 
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GRB peak spectral energy and the isotropic-equivalent/jet 
collimated energy.] 

For objects distributing from z ~ 0.1 to z ~ 6, it is hard 
to believe that they do not evolve. There are cumulative ev- 
idences suggesting that long-duration GRBs prefer to occur 
in low-metallicity galaxies (Fynbo et al. 2003; Hjorth et al. 
2003; Le Floc'h et al. 2003; Sollerman et al. 2005; Fruchter 
et al. 2006; Stanek et al. 2006). With a sample of five nearby 
GRBs, Stanek et al. (2006) have found that the isotropic en- 
ergy of GRBs is anti-correlated with the metallicity in the 
host galaxy (see, however, Wolf & Podsiadlowski 2007). It is 
well known that metallicities evolve strongly with the cos- 
mological redshift (Kewley & Kobulnicky 2005; Savaglio et 
al. 2005). Hence, the evolution of GRBs with the redshift is 
naturally expected (see, e.g., Langer & Norman 2006). 

In this Letter, we use the Amati relation as an example 
to test the cosmic evolution of GRBs. The Amati relation 
is a correlation between the isotropic-equivalent energy of 
long-duration GRBs and the peak energy of their integrated 
spectra in the GRB frame (Amati et al. 2002) 



log E iBO = a + b log £ peak 



(1) 

is defined in the 1— 



The isotropic-equivalent energy E- la 
10000 keV band in the GRB frame. 

With a sample of 41 long GRBs with firmly determined 
redshifts and peak spectral energy, Amati (2006) has ob- 
tained that a = —3.35 and b = 1.75 with the least squares 
method (-E pea k in keV and Ei so in 10 52 erg); and a = -4.04 
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and b = 2.04 with the maximum likelihood method with an 
intrinsic dispersion in the relation (1) being included. Long 
GRBs detected by Swift and having measured redshifts and 
Bpoak are found to be consistent with the Amati relation 
(Amati 2007). 

The difference in the values of the parameters obtained 
with the two methods can be explained as follows. The maxi- 
mum likelihood method directly probes the intrinsic relation 
between the two variables, x = log _Ei so and y = log i? pea k 
(D'Agostini 2005). However, roughly speaking, the least 
squares method estimates the average value of x at a given 
y, (x) = a' + b'y. Teerikorpi (1984) has shown that, when 
x has a Gaussian distribution with a dispersion a x , and the 
relation x = a + by has an intrinsic dispersion a\ in y, b' is 
related to b by 



b' = b 1 + 



(2) 



Amati (2006) has found that cr x « 0.9, b w 2.04, and a\ « 
0.15. Then by equation (2) we have b' w 1.83, which is close 
to the value of 1.75 obtained by the least squares method. 

To test if the Amati relation varies with the cosmolog- 
ical redshift, in this Letter we separate a sample of 48 long 
GRBs [consisting of the long 41 GRBs from Amati (2006) 
and seven additional Swift long GRBs from Amati (2007)] 
into four groups by the GRB redshift. That is, we sort the 
GRBs by their redshifts, and divide them into four groups 
with redshifts distributing from low values to high values. 
Each group contains 12 GRBs (for details see Section 2). We 
then fit each group by equation (1) and calculate the mean 
redshift, and check if the values of a and b evolve with the 
redshift. 

As we will see that, the values of a and b strongly vary 
with the redshift. The variation is not likely to arise from 
the selection effect and hence may indicate that GRBs evolve 
strongly with the cosmological redshift. 

Throughout the Letter, we follow Amati (2006) to adopt 
a cosmology with fi m = 0.3, Oa = 0.7, and Ho = 70 km s _1 
Mpc" 1 . 



2 VARIATION OF THE AMATI RELATION 
WITH THE COSMOLOGICAL REDSHIFT 

To test if the Amati relation (1) evolves with the redshift, 
we separate a sample of 48 long GRBs into four groups ac- 
cording to the redshift of the GRBs, then fit each group with 
equation (1). 

The sample contains 41 long GRBs from Amati (2006, 
Table 1), and seven additional Swift long GRBs from Am- 
ati (2007). The additional seven Swift GRBs are 060115, 
060124, 060206, 060418, 060707, 060927, and 061007. Since a 
GRB sample with z < 0.1 is very incomplete, we select only 
GRBs with z > 0.1 and hence GRB 060218 (z = 0.0331) is 
not included. GRB 060614 is also excluded from our sample 
because of the very large uncertainty in its -E pea k (Amati et 
al. 2007). 

A least squares fit to the 48 GRBs as a single sample 
with equation (1) leads to a = —3.42, b = 1.78, with Xr = 
5.9. The Xr is the reduced x 2 , i e., the x 2 of the fit divided 
by the degree of freedom. A maximum likelihood fit, which 
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Figure 1. Least squares fit to each of the four groups of GRBs 
(data points with error bars; see text) by equation (1) (solid line). 
The two dashed lines mark the l-cr deviation of the fit. 



includes an intrinsic dispersion cri in log -E pca k in the relation 
(1), leads to a = -4.08, b = 2.04, and o\ = 0.14. These 
results are consistent with that obtained with 41 GRBs by 
Amati (2006). 

The redshift of the 48 GRBs spans a range of 0.17-5.6. 
GRB 030329 has the minimum redshift (z = 0.17). GRB 
060927 has the maximum redshift (z — 5.6). The mean red- 
shift is (z) = 1.685. We separate the 48 GRBs into four 
groups with redshifts from low to high, each group contains 
12 GRBs: 

Group A— 12 GRBs, 0.1 < z < 0.84, (z) = 0.56; 
Group B— 12 GRBs, 0.84 ^ z < 1.3, (z) = 1.02; 
Group C— 12 GRBs, 1.3 «S z < 2.3, (z) = 1.76; 
Group D— 12 GRBs, 2.3 < z < 5.6, (z) = 3.40. 

The least squares fit to each group of GRBs by equation 
(1), taking into account the errors in both i5 pca k and -Ei SO , is 
shown in Fig. 1. Immediately one can see that, except Group 
B, the Xr f° r each group is smaller than that obtained by 
fitting the whole sample of GRBs. This fact indicates that 
treating the GRBs at different redshifts as a single sample 
may increase the data dispersion (see Fig. 3 below). 

We find that the values of a and b vary with the mean 
redshift of the GRBs monotonically. 1 In Fig. 2 we plot a and 
b against (z). Clearly, a and b are correlated/anti-correlated 
with (z). The Pearson linear correlation coefficient between 
a and (z) is r(o, (z)) = 0.975, corresponding to a probability 
P = 0.025 for a zero correlation. The correlation coefficient 
between b and (z) is r(b, (z)) — —0.960, corresponding to a 
probability P — 0.040 for a zero correlation. 

A least squares linear fit to a-(z) (the solid line in the 



1 Generally, the variations of a and 
e.g., Li & Paczyriski 2006 



are not independent, see 
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Figure 2. The fitted values of a and b against the mean redshift 
of GRBs. Each data point with error bars represents a group of 
GRBs (A, B, C, and D). The solid line is a least squares linear fit 
to a-(z) and b-{z). 



upper panel of Fig. 2) leads to 

a = -4.58(±0.36) + 0.43(±0.17) z , (3) 

with Xr = 0.13. A least squares linear fit to b-{z) (the solid 
line in the lower panel of Fig. 2) leads to 

b = 2.32(±0.15) - 0.207(±0.066) z , (4) 

with Xr = 0.31. 

The results indicate that a and b strongly evolve with 
the cosmological redshift. 

In Fig. 3 we plot the deviation of fit (s; see Bevington 
& Robinson 1992, Li & Paczyhski 2006) against the mean 
redshift of GRBs. There is not a clear trend for s to vary 
with (z). But it appears that the deviation of fit of each 
group is smaller than that of the whole sample. 



3 IS THE VARIATION CAUSED BY THE 
SELECTION EFFECT? 

To check if the variation of a and b with the cosmological 
redshift is caused by the selection effect, we use Monte-Carlo 
simulations to generate a sample of GRBs according to a pre- 
assumcd Amati relation (1) and with a limit in the observed 
GRB fluence. Then, we divide the sample into four groups 
by the GRB redshift and fit each group by equation (1), just 
as we did in Section 2. 

The lower limit in the bolometric fluence, Fboi.nm, leads 
to a lower limit in the isotropic-equivalent energy of a de- 
tectable burst at redshift z 



E isoAim = AnD 2 coln (l + z)Fi 



bol,lim i 



(5) 



where D com is the comoving distance to the burst. 

In Fig. 4 upper panel, we plot the isotropic energy of the 
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Figure 3. The deviation of fit. Each circle corresponds to a group 
of GRBs. The star represents the result obtained by fitting the 
whole sample (48 GRBs), which is s = 0.16. 



48 GRBs in the sample of Amati (2006, 2007) against their 
redshifts. The isotropic energy is clearly correlated with the 
redshift, with a Pearson linear correlation coefficient r = 
0.437 and a probability P = 0.0019 for a zero correlation. 
The dashed line in the figure is the limit given by equation 
(5) with fboi.iim = 1.2 x 10~ 6 erg cm~ 2 , which reasonably 
represents the selection effect. 

The distribution of the redshifts of the GRBs in the 
sample is plotted in Fig. 4 lower panel. It can be fitted by 
a log-normal distribution, with a mean fi — 0.151 and a 
dispersion a = 0.332 in log z. The Xr of the fit is 0.25. Then, 
the frequency distribution in log z is 



A (log 3) = 



27TO- 



■ exp 



(log z- iif 



2a 2 



(6) 



whose integration over logz (from — oo to oo) is unity. 

The distribution of the isotropic-equivalent energy is 
also described by a log-normal distribution, with a mean 
= 1.09 and a dispersion = 0.85 in log£i so (-Eiso in 10 52 erg). 

Define x = log_Ei so and y = log Speak, where E[ so is 
in 10 erg, and -E pca k is in keV. Assuming that the Amati 
relation is valid and independent of the redshift, and for a 
given x we have y — mx + p with an intrinsic dispersion <j\ 
in y. Then, for a given x, the Gaussian distribution of y is 
given by 



Mv) 



l 



271-CTi 



• exp 



(y-mx-pf 



2a? 



(7) 



By our maximum likelihood fit results in Section 2, we take 
m = 0.49, p = 2.00, and o\ = 0.14. 

The Monte-Carlo simulation is done as follows. First, we 
generate N redshifts with the distribution in equation (6). 
Then, at each redshift, we generate an isotropic-equivalent 
energy with a log- normal distribution (mean of log Ei so — 
1.09, dispersion = 0.85) and satisfying i5i so > Siso.iim- Fi- 
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Figure 4. Upper panel: The isotropic-cquivalent energy versus 
the redshift for the 48 GRBs in the sample. The dashed line is 
the limit given by equation (5) with -Fboljim = 1-2 X 10~ 6 erg 
cm~ 2 . Lower panel: The redshift distribution for the GRBs (his- 
togram). The vertical error bar represents the Poisson fluctuation. 
The dashed curve is a fit to the AT-log z relation by a Gaussian 
function. The bin size in log z is 0.1778. 

nally, for any pair of (z,E iBO ), we generate a peak energy 
Speak according to the distribution in equation (7). Then 
we have a sample of N GRBs, each GRB has a redshift, 
a peak spectral energy, and an isotropic-equivalent energy. 
These GRBs satisfy the distributions described above, and 
the select condition defined by equation (5) (with F bo ui m = 
1.2 x 1CT 6 erg cm -2 ). 

With the above approach, we generated N — 4000 
GRBs. We divided them into four groups by redshift, and 
each group contains 1000 GRBs. Then, we fitted each group 
of GRBs by equation (1) and got the values of a and 6, and 
checked the evolution of a and b with the mean redshift (z). 
We repeated the process 10000 times, each time with a dif- 
ferent sample of 4000 GRBs. We found that a and b indeed 
varied with (z) 

a = -3.57 + 0.105 z , b = 1.84 - 0.0347 z . (8) 

This variation was caused by the selection effect, i.e. the 
limit in equation (5). If we turned off the limit, we found 
that a and b did not evolve with (z). However, comparing 
equation (8) to equations (3) and (4), we found that the 
selection effect is not likely the cause for the evolution in 
equations (3) and (4), since the a and 6 in equation (8) evolve 
too slowly with z. Even if we increased the value of F bo i,ii m 
to 10~ 5 erg cm~ 2 , we got da/dz — 0.22 and db/dz = —0.065, 
whose values are still too small to explain the evolution in 
equations (3) and (4). 

Of course, equation (8) only describes the average evolu- 
tion of a and b for the 10000 runs. For each run, the evolution 
may deviate from equation (8) . To obtain the chance proba- 
bility for the evolution in equations (3) and (4) to arise from 
the selection effect, we used the Monte-Carlo simulation de- 
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Figure 5. The slopes da/dz and db/dz obtained from Montc- 
Carlo simulations (points; 500 runs of 48 GRBs). The cross marks 
the average values of da/dz = 0.105 and db/dz = —0.0347, ob- 
tained with 10000 runs of 4000 GRBs (see text). The star is the 
value obtained with the 48 observed GRBs (eqs. 3 and 4). 

scribed above to generate N — 48 GRBs and repeated the 
process 500 times. For each 48 GRBs obtained in each run, 
we separate them into four groups and calculate da/dz and 
db/dz just as we did to the 10000 runs of 4000 GRBs. The 
results are shown in Fig. 5. 

Based on our simulations (500 runs of 48 GRBs), we 
found that the probability for getting a pair of (da / dz ,db / dz) 
with da/dz > 0.43 and db/dz < -0.207 is 0.04. Hence, we 
have only ~ 4 percent of chance that the variation presented 
in Section 2 is caused by the selection effect. 



4 CONCLUSIONS 

If GRBs do not evolve with the redshift and the selection 
effects are not important, we would expect that the Amati 
relation does not change with the redshift. Hence, from the 
variation of the Amati relation with the redshift we may get 
some clues on the cosmic evolution of GRBs. 

By dividing the 48 GRBs in Amati (2006, 2007) into 
four groups by their redshifts and fitting each group sepa- 
rately, we have found that the isotropic-equivalent energy 
and the peak spectral energy of GRBs remain being cor- 
related in each group, even with a smaller dispersion than 
that for the whole sample. However, the parameters a and 
b in the Amati relation (1) evolve strongly with the redshift 
(eqs. 3 and 4). 

Although the selection effect arising from the limit in 
the GRB fiuence may cause a similar variation of a and b 
(eq. 8), generally the variation is too slow to explain what 
we have found for the observed GRBs (eqs. 3 and 4). With 
Monte-Carlo simulations we have shown that there is only 
~ 4 percent of chance that the observed variation is caused 
by the selection effect. Hence, the variation of the Amati 
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relation with the redshift that we have discovered may reflect 
the cosmic evolution of GRBs and indicates that GRBs are 
not standard candles. 

Our results are limited by the small number of GRBs in 
the sample: we have 48 GRBs in total, and only 12 GRBs in 
each group. To get a more reliable conclusion, the number 
of GRBs with well determined redshifts and spectra need be 
significantly expanded. Since the launch of Swift, the frac- 
tion of GRBs with measured redshifts has increased rapidly. 
However, unfortunately, due to the narrow energy range 
of the Burst Alert Telescope (BAT) on Swift, the fraction 
of bursts that have accurately determined peak/isotropic 
energy has not increased proportionally. The Gamma-ray 
Large Area Space Telescope (GLAST) scheduled for launch 
in late 2007 will provide us with more promise for this pur- 
pose (Omodei 2006). 

We must also stress that our treatment on the selection 
effect has been greatly simplified. The GRBs in the sample 
were detected and measured by different instruments, hence 
the selection effect is much more complicated. A more careful 
consideration of the various selection biases is required to 
determine if the observed evolution of the Amati relation 
reflects the cosmic evolution of GRBs. 

No matter what the conclusion will be (the variation of 
parameters is caused by the GRB evolution effect or by the 
selection effect), our results suggest that it is a great risk to 
use GRBs with redshifts spanning a large range as a single 
sample to draw physics by statistically analyzing the cor- 
relations among observables. Although we have only tested 
the Amati relation, it would be surprising if any of the other 
relations (e.g., the Ghirlanda relation; Ghirlanda, Ghisellini 
& Lazzati 2004a) does not change with the redshift. 
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